Theory of huge tunnehng magnetoresistance in graphene 
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We investigate theoretically the spin-independent tunneling magnetoresistance effect in a graphene 
monolayer modulated by two parallel ferromagnets deposited on a dielectric layer. For the parallel 
magnetization configuration, Klein tunneling can be observed in the transmission spectrum, but at 
specific oblique incident angles. For the antiparallel magnetization configuration, the transmission 
can be blocked by the magnetic-electric barrier provided by the ferromagnets. Such a transmission 
discrepancy results in a tremendous magnetoresistance ratio and can be tuned by the inclusion of 
an electric barrier. 



Recent experiments have demonstrated the stability 
of graphene (a single atomic layer of graphite) and the 
feasibility of controlling its electrical properties by local 
gate voltages ^i^i^i^i^ opening a promising way to ex- 
plore carbon-based nanoelectronics. In graphene, the en- 
ergy spectrum of carriers consists of two valleys labeled 
by two inequivalent points (referred to as K and K') at 
the edges of the hexagonal Brillouin zone. In each val- 
ley, the energy dispersion relation is approximately linear 
near the points where the electron and hole bands touch. 
Such a peculiar band structure results in many interest- 
ing phenomena, including the half-integer quantum Hall 
effecti^'-^^ and minimum conductivity. Further, Dirac- 
like fermions in graphene can transmit through high and 
wide electrostatic barriers almost perfectly, in particu- 
lar for normal incidence j^iSii^ Such a phenomenon, known 
as Klein tunneling, leads to a poor rectification effect in 
graphene p-n junctions^ and thus may limit the perfor- 
mance of graphene-based electronic devices. 

Very recently, inhomogeneous magnetic fields on the 
nanometer scale have been suggested to confine mass- 
less two-dimensional (2D) Dirac electrons, providing 
another clue to the manipulation of electrons in graphene. 
For conventional semiconductor two-dimensional electron 
gas systems, the patterned local magnetic fields define 
various magnetic nanostructures ranging from magnetic 
barriers and well^ to magnetic dots and antidots.^^ 
A great deal of experimental and theoretical works 
have been devoted to understand physical properties of 
Schrodinger fermions in these systems. The effects of 
nonuniform magnetic field modulations on 2D Dirac- 
Weyl fermions, however, has not been investigated as 
thoroughly, especially for the Klein tunneling under inho- 
mogeneous magnetic field. In this work we explore ballis- 
tic transport features of graphene under the modulations 
of both local magnetic fields and local electrostatic bar- 
riers generated by two parallel ferromagnetic stripes. A 
remarkable tunneling magnetoresistance (TMR) effect is 
predicted and its physical mechanism is explained. 

The system under consideration is a single-layer 
graphene sheet covered by a thin dielectric layer, as 
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FIG. 1: (Color online) (a) Schematic illustration of the con- 
sidered two-dimensional electron system modulated by two 
FM stripes deposited on top of the graphene plane. Each FM 
stripe has a rectangular cross section and a magnetization di- 
rected along the current direction (the x axis) . The gate volt- 
age Vg applied on both FM stripes provides an electrostatic 
double barrier in the underneath graphene plane, (b) Simpli- 
fied profiles of the magnetic barrier for the P alignment (spike- 
like lines), the corresponding vector potential Ay[x) (solid 
line), and the electrostatic potential U{x) (dashed line), (c) 
The same as in (b) but for the AP alignment. 



sketched in Fig. [TJa). Two parallel ferromagnetic metal 
(FM) stripes are deposited on top of the insulating layer 
to influence locally the motion of Dirac electrons in the 
graphene (x, y) plane. Both FM stripes have a width 
d and a magnetization in parallel or in antiparallel to 
the current direction (the x axis). Their fringe fields 
thus provide a perpendicular magnetic modulation 5^, 
which is assumed to be homogeneous in the y direction 
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and only varies along the x axis. A suitable external 
in-plane magnetic field can change the relative orienta- 
tion of the two magnetizations which are antiparallel at 
zero field. At the limit of a small distance between the 
graphene plane and the ferromagnets, the magnetic bar- 
rier can be approximated by several delta functions, i.e., 
B,{x) = BIbo{[S{x + L/2)-6{x + D/2)]^j[5{x-D/2)- 
d{x — L/2)]}. Here, B gives the strength of the local 
magnetic field, Ibq = y^h/eBo is the magnetic length for 
an estimated magnetic field 7 represents the mag- 
netization configuration [±1 or parallel (P) /antiparallel 
(AP)], D is the distance between the two FM stripes, and 
L = 2d-\- D is the total length of the structure along the 
transport direction. The model magnetic field configu- 
rations for 7 = ±1 are depicted in Figs. [Ifb) and[ll^c), 
respectively. Further, when a negative gate voltage is 
applied to both FM stripes, a tunable electrostatic dou- 
ble barrier potential U{x) arises in the graphene layer. 
A square shape with height Uq can be taken for the 
electric potential created by either gate. Accordingly, 
the simplified electrostatic barrier has the form, U {x) = 
Uo[e{x + L/2)e{-D/2 - x) + e{x - D/2)Q{L/2 - x)], 
where 0(x) is the Heaviside step function. 

For such a system, the low-energy excitations in the 
vicinity of the K point can be described by the following 
Dirac equation 
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FIG. 2: (Color online) Transmission as a function of incident 
angle for electrons traversing the considered structure (de- 
picted in Fig. [T]) with (a)-(d) parallel or (e)-(f) antiparallel 
magnetization configuration. Device parameters used in the 
calculations are d — D — 1 and B — 2. The incident energy 
is fixed at ^ = 5. Note that all curves in (f) are scaled by a 
factor 200. 



where v^, ~ 0.86 x 10^ m/s is the Fermi velocity of the 
system, (Jx^ (Jy^ and cr^ are three isospin Pauli matrices, 
P = iPx-tPy) is the electron momentum, A is the vec- 
tor potential which in the Landau gauge has the form 
A =(0, A^(x), 0), and ctq is the 2x2 unit matrix. Since 
the Dirac Hamiltonian of graphene is valley degenerate, 
it is enough to consider the K point. For convenience 
we express all quantities in dimensionless units by means 
of two characteristic parameters, i.e., the magnetic length 
Ibq and the energy Eq = hv^/lBo- For a realistic value 
5o = 0.1 T, we have /^o = 811 A and Eq 7.0 meV. 

Since the system is homogeneous along the y direction, 
the transverse wave vector ky is conserved. At each re- 
gion with a constant vector potential Ay and electrostatic 
potential /7, the solution of Eq. ([1]) for a given incident 
energy E can be written as 
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Here q = ky-\- Ay^ and kx is the longitudinal wave vector 
satisfying 



kl + {ky+Ayf = {E-Uf. 



(3) 



The sign of kx is chosen in such a way that the corre- 
sponding eigenstate is either propagating or evanescent in 
the forward direction. The coefficients C+ and C- are de- 
termined by the requirement of wave function continuity 
and the scattering boundary conditions. The scattering 



matrix method^^ is adopted to obtain these coefficients 
and the transmission probability T = Tp^j^p{E^ ky) for 
the P/AP configuration. The latter depends on the in- 
cident energy E and the transverse wave vector ky. The 
ballistic conductance at zero temperature is calculated 
from 
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where Ly ^ L is the sample size along the y direction, 
6 is the incident angle relative to the x direction, and 
Go = 2e^EFLy/ {t^K) is taken as the conductance unit. 

The proposed device relies on the interplay between 
the Klein tunneling and the wave vector filtering pro- 
vided by local magnetic fields. To obtain a quantitative 
understanding of this interplay. Fig. [2] plots the transmis- 
sion probability calculated as a function of the incident 
angle for both the P and AP magnetization configura- 
tions. In our calculations the structure parameters of the 
magnetic barrier are sei si d = D = 1 and B = 2. The 
incident energy is fixed at £^ = 5 and the electric barrier 
height Uq is taken as 2,3, 4, 5, 6, 7, 8 for different curves. 

For the magnetic barrier with P alignment, the 
transmission spectrum demonstrates obvious angular 
anisotropy [see Figs. [2](a)-[2l^d)]. The reffection at nor- 
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mal incidence is finite and is almost complete at suitable 
electric barrier heights. Instead, perfect transmission ap- 
pears at some oblique incidences. For example, in the 
special case E = Uq^ the transmission peak with a fi- 
nite width locates at ky = —Ay [see Fig. [2fd)]. In com- 
parison with the case of pure electric barriers,^ we can 
see that the magnetic barrier changes the incident direc- 
tion at which the Klein tunneling occurs. The transmis- 
sion is remarkable in a wide region of negative and is 
blocked by the magnetic barrier when the incident an- 
gle exceeds a critical value 0-^{Uo) or is below another 
critical value 6-{Uo). This can be understood as follows. 
From Eq. ([3|) we know that evanescent states appear in 
the magnetic barrier regions when the magnetic vector 
potential (here Ay = B) and electrostatic barrier sat- 
isfy \ky ^ B\ > \E — Uo\. The transmission is generally 
weak as the decaying length of the evanescent states is 
shorter than the barrier width. In the transmission for- 
bidden region, there may exist one or two line-shaped 
peaks with unity values, as a result of resonant tunnel- 
ing through the symmetric double barrier structure. The 
applied electric barrier significantly alters the positions 
of the transmission peaks. We can also observe a large 
difference between the transmission curves for the bar- 
rier height Uq < E and Uq = 2E — Uq. Such a difference 
arises from the fact that the carrier states for the two 
cases are not completely complementary. 

We next examine the transmission characteristics for 
the AP alignment, which is shown in Figs.[2](e) and[2](f). 
In this configuration the magnetic vector potential is an- 
tisymmetric about the central line x = [see Fig. [U^c)]. 
The Dirac Hamiltonian possesses a symmetry associated 
with the operation TRx&y^ where T is the time rever- 
sal operator and Rx is the reflection operator x —x. 
This symmetry implies the invariance of the transmis- 
sion with respect to the replacement ky —ky^ as seen 
in Figs. [2je) and[2](f). For large \E — Uo\ the transmis- 
sion decays monotonically as the incident angle increases 
from zero [see Fig. El^e)]. Since the carrier states in the 
two magnetic barriers are not identical, perfect transmis- 
sion can not be achieved (except for the case of normal in- 
cidence). Note that for the AP configuration and a given 
wave vector /c^ ^ 0, the presence of evanescent states in 
the first magnetic barrier only requires ky > \E — Uo\—B. 
When \E — Uo\ < B this condition is met for all in- 
cident directions and the transmission can be strongly 
suppressed, as shown in Fig.lSff). Within this parameter 
regime, the transmission exhibits a nonmonotonic varia- 
tion with the positive incident angle. Furthermore, the 
maximal transmission for the AP alignment can be 2 or- 
ders of magnitude lower than that for the P alignment 
[Figs.EJc) andEJd)]. 

As demonstrated above, the transmission features for 
the P and AP configurations are quite distinct. Such 
a difference is also exhibited in the measurable quan- 
tity, the conductance G. In Fig. [3] the conductance is 
plotted as a function of the Fermi energy for several 
heights of the electric barrier. Resonant peaks can be 
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FIG. 3: (Color online) Conductance as a function of Fermi 
energy for electrons traversing the considered structure with a 
parallel (solid line) or antiparallel (dashed line) magnetization 
configuration. Device parameters used in the calculations are 
d = D = 1 and B = 2. 
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FIG. 4: (Color online) MR ratio as a function of [(a) and 
(c)] Fermi energy or [(b) and (d)] electric barrier height for 
electrons traversing the considered magnetic-electric barrier 
structure. In (a) and (b) simplified magnetic field profiles are 
utilized and the device parameters used in the calculations are 
d = D = 1 and B = 2. In (c) and (d) realistic magnetic field 
profiles are taken. In the calculations we assume that both 
ferromagnetic stripes have a rectangular cross section of width 
d = 1 and height dz =0.6 and magnetization iiqMx = 1.8 T 
(for cobalt material), placed at a distance of zq — 0.2 on top 
of the graphene plane. Their distance is D = 1. 



observed in the conductance spectrum for both P and 
AP alignments. For the P alignment, the conductance is 
finite (larger than 0.1 in most cases in Fig. [3|). For the 
AP alignment, the conductance is almost zero within a 
broad energy interval [covering {Uq — B^Uq ^ B)] except 
for several sharp conductance peaks. In this energy re- 
gion Gap is depleted by the magnetic barrier whereas 
Gp is finite. Away from this transmission-blocking re- 
gion. Gap essentially increases with the Fermi energy 
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and is primarily contributed by the propagating modes. 
The normahzed difference between Gp and Gap^ i-e., 
the TMR ratio MR = {Gp- Gap) /Gap, is presented in 
Fig. [41(a). In the absence of the electric barrier, high val- 
ues of MR are located in the low Fermi energy region, as 
a result of the strong suppression of transmission in the 
AP alignment. The inclusion of an electric barrier shifts 
the transmission-blocking region and, thus, can be used 
to adjust the MR ratio. The latter is obviously reflected 
in Fig.H^b). 

In the above analysis, we take simplified magnetic field 
profiles to illustrate the operating principles of the pro- 
posed device. In realistic cases the modulated magnetic 
field Bz{x) has the smoothing variations on the scale of 
graphene lattice spacing (a = 0.246 nm). When both FM 
stripes have the same rectangular cross section and mag- 
netization along the x-direction, the generated magnetic 
field profiles for the P and AP alignments can be ob- 
tained analytically.^^ For the parameters given in the fig- 
ure caption the calculated MR ratio is shown in Figs. [31(c) 
and [11(d). The calculation shows that the conductance of 
the device has a variation similar to that in Figs. [41(a) 
and [11(b). The TMR ratio remains large and exhibits 
rich variations as the electric barrier height increases. 
Since ferromagnetic elements with a submicron scale have 
been successfully fabricated on top of a two-dimensional 
electron system^^ and dielectric layers on monolayered 
graphene have been realized recently, "^^^ our considered 
structure is realizable with current technology. 



In summary, we have investigated the transport fea- 
tures of a graphene monolayer under the modulation of 
both a magnetic double barrier and an electric barrier, 
where the magnetic double barrier is provided by deposit- 
ing two parallel ferromagnetic stripes with magnetiza- 
tions along the current direction. The results indicate 
that for the AP magnetization configuration the trans- 
mission of electrons in graphene can be drastically sup- 
pressed for all incident angles. When in the P alignment 
the Klein tunneling can be generally observed at specific 
oblique incident directions rather than the normal inci- 
dence. The difference of wave-vector-dependent trans- 
mission for two magnetization configurations (P/AP) 
leads to a large TMR ratio, which can be further adjusted 
by the electric barrier. Note that different thin dielec- 
tric layers atop graphene sheets have been fabricated and 
then the top gates can be formed by means of standard 
e-beam lithography.^^ ^ The deposition of ferromagnetic 
materials on insulating layers has been widely adopted to 
create local magnetic field modulations of the underlying 
2D semiconducting systems.—'— Thus our proposed de- 
vice is within the realizable scope of current technological 
advances. 
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